APPROXIMATION OF QUANTUM GRAPH VERTEX COUPLINGS 
BY SCALED SCHRODINGER OPERATORS ON THIN BRANCHED 

MANIFOLDS 



PAVEL EXNER AND OLAF POST 



OO 

o 



Abstract. We discuss approximations of vertex couplings of quantum graphs using 
families of thin branched manifolds. We show that if a Neumann type Laplacian on such 
manifolds is amended by suitable potentials, the resulting Schrodinger operators can 
, approximate non-trivial vertex couplings. The latter include not only the (5-couplings 

' but also those with wavefunctions discontinuous at the vertex. We work out the example 

^ , of the symmetric (5'-couplings and conjecture that the same method can be applied to 

Q ' all couplings invariant with respect to the time reversal. 

(N 
(N 

1. Introduction 

^ • The quantum graph models represent a simple and versatile tool to study numerous 

physical phenomena. The current state of art in this field is described in the recent 



proceedings volume |EKK+08j to which we refer for an extensive bibliography. 

One of the big questions in this area is the physical meaning of quantum graph vertex 
coupling. The general requirement of self-adjointness admits boundary conditions con- 
taining a number of parameters, and one would like to understand how to choose these 
when a quantum graph model is applied to a specific physical situation. One natural 
' idea is to approximate the graph in question by a family of "fat graphs", i.e. tube-like 

o ■ 
m 



O . manifolds built around the graph "skeleton" , equipped with a suitable second-order dif- 
^ ferential operator. Such systems have no ad hoc parameters and one can try to find what 
boundary condition arise when the manifold is squeezed to the graph. 

The question is by no means easy and the answer depends on the type of the operator 
^ ■ chosen. If it is the Laplacian with Dirichlet boundary conditions one has to employ 
an energy renormalisation because the spectral threshold given by the lowest transverse 
eigenvalue blows up to infinity as the tube diameter tends to zero. If one chooses the 
^ ' reference point between the thresholds, the limiting boundary conditions are determined 
■ by the scattering on the respective "fat star" manifold |MV07j . If, on the other hand, 
the threshold energy is subtracted, the limit gives generically a decoupled graph, i.e. the 
family of edges with Dirichlet conditions at their endpoints |P05t IMV071 IDTOGj . One 
can nevertheless get a non-trivial coupling in the limit if the tube network exhibits a 
threshold resonance |G08lRCF07j . and moreover, using a more involved limiting process 
one can get also boundary conditions with richer spectral properties |CE07j . 

The case when the fat graph supports a Laplacian of Neumann type is better un- 
derstood and the limit of all types of spectra as well as of resonances has been worked 
out |FW93l [RSOTl IKuZOll [EFOSl lEPOTl iGOSl [EP08] . Moreover, convergence of resol- 
vents etc. has been shown in |SaiOOl \POQ\ IEP07] . Of course, no energy renormalisation 
is needed in this case. On the other hand, the limit yields only the simplest boundary 
conditions called free or Kirchhoff. 

The aim of this paper is to show that one can do better in the Neumann case if the 
Laplacian is replaced by suitable families of Schrodinger operators with properly scaled 
potentials. Such approximations have been shown to work on graphs themselves |E97| 
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lENZOl] , here we are going to "lift" them to the tube-hke manifold£|. First we will show 
that using potentials supported by the vertex regions of the manifold with the "natural" 
scaling, as where e is the tube radius parameter, we can get the so-called 5-coupling, 
the one-parameter family with the wavefunctions continuous everywhere, including at 
the vertex. This shows, in particular, that one cannot achieve such an approximation in 
a purely geometric way, with a curvature-induced potential of the type |DEK01] , because 
the latter scales typically as e^"^. As main result in this case, we show the convergence 
of the spectra and the resolvents (cf. Theorems l3.3ti3Tl) . 

On the other hand, the 5-coupling is only a small part in the set of all admissible 
couplings; in a vertex joining n edges the boundary conditions contain ra^ parameters. 
Here we use the seminal idea of Cheon and Shigehara |CS98j applied to the graph case 
in |CE04j and generalised in |ET06t IET07] . For simplicity we will work out in this paper 
the example of the so-called symmetric ^'-coupling, in short Si,, a one-parameter family 
which is a counterpart of 6, by using the result of |CE04] and "lifting" it to the manifold. 
We show that such a coupling is approximated with a potential in the vertex region 
together with potentials at the edges with compact supports approaching the vertex, 
all properly scaled, cf. Theorem 14. 7[ The speed with which the potentials are "coming 
together" must be slower than the squeezing; the rate between the two we obtain is 
surely not optimal. 

We have no doubts that in the same way one can lift to the manifolds the more general 
limiting procedure devised in |ET07j which gives rise to a ("^^)-parameter family of 
boundary conditions, namely those which are invariant with respect to the time reversal. 
Such an extension would be technically demanding, however, and in order not to burden 
this paper with a complicated notation and voluminous estimations we postpone it to a 
later work. 

Let us survey the contents of the paper. In the next section we define the graph and 
manifold models and provide necessary estimates. In Section 12] we prove the convergence 
in the 5-coupling case. For simplicity we restrict ourselves to a star-shaped graph with 
a single vertex; the approximation bears a local character and extends easily to more 
complex graphs. Finally in Section 4 we extend the result to the ^g-coupling case and 
comment on the applicability of the method to more general couplings. 



2. The graph and manifold models 

2.1. Graph model. Let us start with a simple example of a star-shaped metric graph 
G = ly having only one vertex v and degf adjacent edges E = Ey of length 4 ^ (0, oo], 
so we can think of E = {1, . . . , degf }. We identify the (metric) edge e with the interval 
le '■= (0,£e) oriented in such a way that corresponds to the vertex v. Moreover, the 
metric graph G = is given by the abstract space ly := 1^/ ~ where U denotes the 
disjoint union, and where the equivalence relation ~ identifies the points G /g with the 
vertex v. The basic Hilbert space is \-2{G) := ^^^^ '-2(-^e) with norm given by 

ii/ir = ii/iiG = E ri/(^)i'd^- 

The decoupled Sobolev space of order k is defined as 

HLx(G) :=0H'=(/e) 



"'^ Another approach to approximation of nontrivial vertex conditions was proposed recently in |Pa07b| 
IPa07a) 
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together with its natural norm. Let p = {pe}e be a vector consisting of the weights 
Pe > for e G E. The Sobolev space associated to the weight p is given by 

H^(G):={/gHLx(G)|/gC^}, (2.1) 

where / := {/e(0)}e € C'^^^'' is the evaluation vector of / at the vertex v and Cp is the 
complex span of p. We use the notation 

l=f{v)p, i.e., fe = fiv)pe (2.2) 

for all e E E. In particular, if p = (1, . . . , 1), we arrive at the continuous Sobolev space 
H^(G) := Hp(G'). The standard Sobolev trace estimate 

\9m'<a\\gTm + l\\9\\m (2-3) 

for g G H^(0,£) and < a < i ensures that Hp(G) is a closed subspace of y^laaxi^)^ ^.nd 
therefore itself a Hilbert space. A simple consequence is the following claim. 

Lemma 2.1. We have 

\fiv)\'<\pr{a\\f% + l\\f\\l) 
for f G Hp(G) and < a < i_ := miUe^E^e, the minimal length at the vertex v. 

We define various Laplacians on the metric graph via their quadratic forms. Let us 
start with the (weighted) free Laplacian A^. defined via the quadratic form d = Do given 
by 

Hf) ■■= mi = E ll/elll and domO := H^(G) 

e 

for a fixed p (the forms and the corresponding operators should be labelled by the weight 
p, of course, but we drop the index, in particular, because we are most interested in the 
case p = (!,...,!)). Note that is a closed form since the norm associated to the 
quadratic form d is precisely the Sobolev norm given by ||/|Ihi(g) ~ II/'IIg + II/IIg- 

The Laplacian with S-coupling of strength q is defined via the quadratic form f) = f)(G,g) 
given by 

m--=\\f'\\l + qiv)\fiv)\' and domf):=H^(G). (2.4) 
The 5-coupling is a "small" perturbation of the free Laplacian, namely we have: 

Lemma 2.2. The form i){G,q) is relatively form-bounded with respect to the free form dc 
with relative bound zero, i.e., for any r] > there exists > such that 



m)-^if)\ = \Qiv)\\fiv)\'<vm + C, 



2 

G- 



Proof. It is again a simple consequence of Lemma 12.11 Since we need the precise be- 
haviour of the constant C,,, we give a short proof here. From Lemma [2.11 we conclude 
that 

\w)-m\<\qmp\-Ham + -mi\ 



for any < a < Set a := min{r7|pp/|g(t>)| , £_} and 

\q{vr \q{v)\ 



then the desired estimate follows. □ 

One can see that the norms associated to f) and d are equivalent and, in particular, 
setting = 1/2 in the above estimate yields we get: 



4 



PAVEL EXNER AND OLAF POST 



Corollary 2.3. The quadratic form f) is closed and obeys the estimate 

()(/)< 2(f)(/) + Cv2||/||^). 

The operator H = H(^G,q) associated to f) acts as {Hf)^ = — on each edge and 
satisfies the conditions 

l^ = l^=:f{v) and J^P^fM = <liv)m (2.5) 

for any pair (ei, 62) of edges meeting at the vertex v. We use the formal notation 

H = H(^G,,) = Aa + qiv)6,; (2.6) 

note that is a non-negative operator by definition. In order to compare the "free" 
quadratic form with the graph norm of H we need the following estimate: 

Lemma 2.4. We have 

II/IIhmg) = Hf) + 11/11^ < 2max{Ci/2, v^}||(i/ - i)/||^ 

for f e domif C domf) = H^(G'). 



Proof. Using the estimate of Corollary 12.31 we obtain 

f (/) + ll/f < 2(f3(/) + (Cv2 + imr) < 2|[)(/) + ll/f I + 2Cv 
Moreover, the first term can be estimated as 

Hf) + ll/f r < 2(f3(/f + ii/r) = 2|f)(/) - iii/irr = 2|(/, (i/ - i)/) p 

<2ii/ini(if-i)/f. 

Finally, we apply the estimate ||/|| < ||(-f^ — i)/|| to obtain the result. □ 

Remark 2.5. Note that we have not said anything about the boundary conditions at the 
free ends of the edges of finite length if there are any. As we employ the Sobolev space 
Hp(G) for the domain, we implicitly introduce Neumann conditions for the operator, 
f'ei^e) = 0. However, one can choose any other condition at the free ends, or to construct 
more complicated graphs by putting the star graphs together. 

2.2. Manifold model of the "fat" graph. Let us now define the other element of 
the approximation we are going to construct. For a given e G (0,eo] we associate a 
(i-dimensional manifold to the graph G in the following way. To the edge e & E and 
the vertex v we ascribe the Riemannian manifolds 

X^,e := h X eY^ and X^^^ := eX^, (2.7) 

respectively, where sY^ is a manifold Ye equipped with metric he^e '■= £^he and eXs^v 
carries the metric ge,v = ^'^Qv here he and are e-independent metrics on Yg and X^, 
respectively. We identify the boundary component deX^^v = sdeX^ of dX^^v = sdX^ with 
dyX^^e = ^d^Xe = {0} X eYe and call the resulting manifold X^. We refer to the unsealed 
manifold (which conventionally means e = 1) as to X. In particular, the manifold X 
consists of the number degf of straight cylinder^ with cross-section Y^ and a vertex 
neighbourhood manifold X^ containing the boundary components Y := IJg Ye- Without 
loss of generality, we may assume that each cross-section Y^ is connected, otherwise 
we replace the edge e by as many edges as is the number of connected components. 
We denote the boundary component of dX^ at the edge e by dgXy and the boundary 
component of dXe at the vertex v by dyXe = {0} x Yg. Note that these two boundaries are 
identified in the entire manifold X. Similarly, we denote by dgX^^^ = edeXy and dyXi^^e = 



^The straightness here refers to the intrinsic geometry only. We do not assume in general that the 



manifolds are embedded, for instance, into a Euclidean space, see also Remark 1 2. 6 
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ed^Xf, their scaled versions. For convenience, we will always use the e-independent 
coordinates (s,?/) G = Jg x Ye and x G X^,, so that the radius- type parameter e only 
enters via the Riemannian metrics. 

The entire manifold X^ may or may not have a boundary c^Xg, depending on whether 
there is at least one finite edge length £e < cxd or one "transverse" manifold has a 
non-empty boundary. In such a situation, we assume that Xg is open in X^ = X^ U SX^. 
A particular case is represented by embedded manifolds which deserve a comment: 

Remark 2.6. Note that the above setting contains the case of the ^-neighbourhood of an 
embedded graph G C M^, but only up to a longitudinal error of order of e. The manifold 
Xe itself does not form an ^-neighbourhood of a metric graph embedded in some ambient 
space, since the vertex neighbourhoods cannot be fixed in the ambient space unless one 
allows slightly shortened edge neighbourhoods. Nevertheless, introducing e-independent 
coordinates also in the longitudinal direction simplifies the comparison of the Laplacian 
on the metric graph and the manifold, and the error made is of order of 0(e), as we will 
see in Lemma [2.71 for a single edge. 

The basic Hilbert space of the manifold model is 

^2{Xe) = 0(L2(/e) ® ^(5^)) © ^2{eX,) (2-8) 
e 

with the norm given by 

where dxe = di/e ds and dx^ denote the Riemannian volume measures associated to the 
(unsealed) manifolds Xe = Je x Yg and X^, respectively. In the last formula we have 
employed the appropriate scaling behaviour, dx^^e = e'^~^dyeds and dx^^v = s'^dxy. 

Denote by H^(Xe) the Sobolev space of order one, the completion of the space of 
smooth functions with compact support under the norm given by ||m||hi(X£) ~ II'^'^IIxe + 
As in the case of the metric graphs, we define the Laplacian A^^^ on X^ via its 
quadratic form 

Mu):=\\du\\l^ = J2^'~' I {\u\s,yr + ^^\dyMl)^yeds + e'-' j | dn|J^ dx. 

(2.9) 

where u' denotes the longitudinal derivative, u' = dsU, and du is the exterior derivative 
of u. As before, the form is closed by definition. Adding a potential, we define the 
Hamiltonian as the operator associated with the form f)^ = i){Xs,Qs) given by 

where the potential has support only in the (unsealed) vertex neighbourhood X. and 

Qe{x) = ^Q{x) (2.10) 

where Q = Qi is a fixed bounded and measurable function on X„. The reason for 
this particular scaling will become clear in the proof of Lemma [3^ Roughly speaking, it 
comes from the fact that volX^^. is of order e"^, whereas the (d— l)-dimensional transverse 
volume vol Fe^e is of order e"^"^. The operators i^e and Ag are associated to forms f)^ and 
Oe, respectively; note that Ag = A^^^ > is the usual (Neumann) Laplacian on X^. As 
usual the Neumann boundary condition occurs only in the operator domain if dXg, 7^ 0. 
We postpone for a moment the check that if^ is relatively form-bounded with respect to 
Aj^ , see Lemma [2.101 below. 
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Let US compare the two cylindrical neighbourhoods, X^^e = / x eYe and X^^e = ^^e, 
on edges of length £ > and £e = (1 — e)i, respectively. The result for the entire space 
Xe then follows by combining the estimates on the edges and the fact that the potential 
is only supported on the vertex neighbourhoods. The verification of the conditions of 
(5-closeness in the next lemma is straightforward; for more details on 5-closeness we refer 
to [P06l App. A] or jP08] . 

Lemma 2.7. Let := L2{X^^e) o-nd M'f. := L2(X£ e)- Moreover, define 

Je-.Ji^e^M {JefKs, y) .= /((I - y) , 

j^: := /((I - e)s,y). 

Then the quadratic forms fe(/) := H/Hx^^ o-nd ^^(m) := with domO^ = H^(X£ 



and domDe = H-'^(Xe^e) o^f^ de-close with 5^ = 2e/{l — ef-l'^; namely, we have J'^J^ = id, 

JeJ^ = id, II Jell < 1,' II Jell < 1 + 4, 

||J:- J:II <4 and |^,(JJ,«)-c),(/,J»| <5,. 
In particular, we get 

m^^^ + - Je(A^^,^ + 1)-V:|| < 24 = 0(e). 

Before we check the closeness assumptions of jPOGj Appendix] in the next section, we 
need some more notation and estimates. The estimates are already provided in |EP05t 
IP06j . but we will also need a precise control of the edge length, when we approximate 
the 5g-coupling by 5-couplings in Section H] below. Therefore, we present short proofs of 
the estimate here. 

We first introduce the following averaging operators 

-f^u:=+ udx^ and j u{s) := i u{s,-)dye 

J X.V J Ye 

for u G L2(Xe), where we use the following symbols 

uAx^:= — [ udxy and / ipdye:= — 7-~~^ f ^^V^ 

X„ void X^ Jx, J Ye Vold-l Ye J Ye 

denoting the normalised integrals. For brevity, we also omit the measure and write 
etc. 

In order to obtain the below Sobolev trace estimate (12.111) . we need a further decom- 
position of the vertex neighbourhood Xy. Recall that Xy has (degf)-many boundary 
components isometric to Yg. We assume that each such boundary component has a collar 
neighbourhood X^^e = (0,-^e) x Ye of length ie- Note that the scaled vertex neighbour- 
hood Xir^y = eXy Is of ordcr e in all directions, so that the scaled collar neighbourhoods 
^e,v,e '■= £^v,e are of length eie- We can always assume that such a decomposition exists, 
by possibly using a different cut of the manifold X into X^ and Xe, the price being an 
additional longitudinal error of order e (see Lemma [2.71) . Similarly as in (12.31) . one can 
get the following Sobolev trace estimates for the scaled manifolds: 

_ 2 

11^11?, V. < ealldwlli +— IImIIy (2.11) 

lhlllx.,.<«lh'lli,. + ^lhllx,e (2.12) 

for < a,a < £e on the vertex and edge neighbourhood, respectively, where u' = dsU 
denotes the longitudinal derivative. The unsealed versions are obtained, of course, by 
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setting e = 1. Moreover, by the Cauchy-Schwarz inequality we get 

vol Fe 1/^^(0)1^ < \\u\\l^x, = \\u\\l^x,. 

In the following lemma we compare the averaging over the boundary of with the 
averaging over the whole space X^: 

Lemma 2.8. For u G H"'^(Xt,), we have 

yo\dX,\^ - ^ u'\ < ^vol9eX,|/g^^^u-f^M|' < («^+^w3)l|d«llL 



eG-B 



/or <a<£_= mine4; where X2iv) denotes the second (i.e., first non-zero) eigenvalue 
of the Neumann Laplacian on ; the latter is defined conventionally as the operator 
associated to the form dy{u) := lldwHl^^ with the domain domO^, := H-'^(X^). 

Proof. Using Cauchy-Schwarz twice and the estimate (12. lip for each edge e and e = 1, 
we obtain 

voldX^\fg^w\^ < ^vol^eX^lfg^^^wj^ < \\u\\lx, < a||dw;||^^ + ^\\w\\l^ (2.13) 

e 

for < a < using the fact that IJe"^f,e C X^. We apply the above estimate to the 
function w = PyU := u — f^u and observe that 



\w\ 



\2 ^ ^ IIJ„..I|2 



^^^W^^Wx. (2-14) 



as one can check using the fact that that dw = du and that is the projection onto 
the orthogonal complement of the first eigenfunction 1^, G L2(X„). □ 

We also need an estimate over the vertex neighbourhood. It will assure that in the 
limit e —>■ 0, no family of normalised eigenfunctions (Me)^ with eigenvalues lying in a 
bounded interval can concentrate on 



Lemma 2.9. We have 



\u\\x < 4:6^ 



+ Cvoi (a + 



2 

II /ii2 II li2 

a\\u \\y H — \\u\\ X 

II ll-^E.B qII ''^£,E 



-A2(f) V aA2(f)' 

for < a,a < = mine4; where Cvoi := Cvoi(f) = volX^,/ voldX^ and X^^e '■= Ue^^e 
denotes the union of all edge neighbourhoods. 

Proof. We start with the estimate 

<2e'^(\\u-{ u\\l +\\-f u\\l] <2e'^(-^\\du\\l +volXJf m'^ 



V 



using Lemma 12.81 and the fact that j^^u is constant. Moreover, the last term can be 
estimated by 



vol9X,|_f^M|' < 2vo\dX,(\f^u-fg^u\^ + \j^g^^ 



u\ 



- + $^vol9,X.|f 



deXv I 



using (12.131) . Since dgX^ is isometric to d^X^, we can estimate the latter sum by 

2 



due to (I2.12P for e = 1, each edge e and < a < Here, Xe '■= Xi^e is the union of 
the unsealed edge neighbourhoods. The desired estimate then follows from the scaling 



PAVEL EXNER AND OLAF POST 



d-1 1 



w\ 



for w = u or w = u' (where 

□ 



behaviour ||dM||Y =e \\du\\y and \\w\\x =e 

II WJ^^^v II II^t; II ll^£,e 

u' = dgU denotes the longitudinal derivative). 

We are now able to prove the relative (form-)boundedness of the Hamiltonian if^ 
with respect to the Laplacian Aj^^ for the indicated class of potentials. It is again 

important here to have a precise control of the constants and Cr, in terms of the 
various parameters of our spaces. This will be of particular importance when we deal 
with the approximation of the 5g-coupling by 5-couplings with shrinking spacing a = 
in Section m below. 

Lemma 2.10. To a given rj G (0, 1) there exists > such that the form i)^ is relatively 
form-bounded with respect to the free form ds with relative bound r] for all e G (0,ej, in 
other words, there exists C,, > such that 



\^e{u) -'Oe{u)\ <ri'i)e{u)+C^r,\\u\ 



whenever Q < £ < Sn, where the constants Etj and are given by 



+ Cvol 



Crj = C',,(||Q||oo,^-) := SCvoillQlloomax 
Note that = 0(i) and C„ = OUZ') as £_ 0. 



4cvoi||g||oo 2 

r] 



(2.15a) 
(2.15b) 



— y\. J u,xxv_j. — 0(^_ ) as £— 

Proof. The potential = £~^Q is by assumption supported on the vertex neighbour- 
hood Xy, therefore we have 



Mf)-Mf)\ < 
<4||g||oc 



WQW 



°o II Il2 
— M 



+ Cvol 



1- + 



|dn"2 



X, 



I II 'Il2 1 



ooCvol I 



U\ 



X, 



E.E 



using Lemma [2^ for < a < £_ and a := Choosing a = min{£_, r7(4cvoi||<5||oo) ^} 
and < e < Erj with as above, we can estimate the quadratic form contributions by 

v{\\M\x,,, + \W\\x,,e) ^ vWMlx,- 

The expression for C„ then follows by evaluating the constant in front of the remaining 



norm. 



□ 



We need to estimate the "free" quadratic form against the form associated with the 
Hamiltonian: 

Corollary 2.11. The quadratic form f)^ is closed. Moreover, setting rj = 1/2, we get the 
estimate ^ 

c),(m) <2([),(m) + Ci/2||m||^J 
which holds provided < £ < £1/2- 

As in Lemma 12. 4^ we can prove the following estimate in order to compare the "free" 
quadratic form with the graph norm of H^: 

Lemma 2.12. We have 

I|w|Ihi(x,) =fe(M) + h||^, <2max{(7i/2,V2}||(i/,-i) 

for u G domi/^ C domf)^ = H^(X£) and Q < e < Eq. 



X, 
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3. Approximation of (5-couplings 

After this preliminaries we can pass to our main problems. The first one concerns 
approximation of a 5-coupling by Schrodinger operators with scaled potentials supported 
by the vertex regions. For the sake of simplicity most part of the discussion will be done 
for the situation with a single vertex as described in Section 12.11 



3.1. Quasi-unitary operators. First we define quasi-unitary operators mapping from 
^ to M' and vice versa, as well as their analogues on the scales of order one, namely 
and Here, 

^ := 12(6-), := H^(G), ^^{X,), := H^(A,). (3.1) 

Moreover, we need a relation between the different constants of the graph and the man- 
ifold model introduced above. Specifically, we set 

Pe := (volrf_iye)'/' and q{v) = I Q^x,. (3.2) 

Let J: — > be given by 

J/:=e-('^-i)/2 0(/^®4e)©O (3.3) 

with respect to the decomposition (12.81) . Here is the normalised eigenfunction of 
associated to the lowest (zero) eigenvalue, i.e. ie{y) = (yold-iY^)'^^'^. In order to relate 
the Sobolev spaces of order one we need a similar map: we define : — > by 

Ji/ := e-('^-i)/2(0(/, ® 4,) © /(^;)1,), (3.4) 

eeE 

where 1^ is the constant function on X„ with value 1. Note that the latter operator is 
well defined: 

(jV)e(0,2/) = e-('^-^)/Ve~Ve(0) = S-^"-'^ f (v) = {j' fUx) 

for any x G X^, due to (13.21) and (12. 2p . i.e., the function J^f matches along the different 
components of the manifold, thus Jf G H^{Xs;). Moreover, f{v) is defined for / G H^(G') 
(see Lemma [2.11) . 

The mapping in the opposite direction, J': Jf — > Jif, is given by the adjoint, J' : = 
J*, which means that 

(J'«)e(5) =£('^-^)/'(4e,ne(5,-))y. =^^'''^/We«(^)- (3-5) 

Furthermore, we define J'^ : J>^^ — > Jf^ by 

{J:\){s) := e^'-^y' u^is, ■))y. + Xe{s)pe - fu{0))] . (3.6) 

Here Xe is a smooth cut-off function such that Xe(0) = 1 and Xe(4) = 0. If we choose 
the function Xe to be piecewise affine linear with Xe(0) = 1, Xe(a-) = and Xe(^e) = 0, 
then ||Xe||/g = ct/S < a and ||Xelll = (^'^^ Moreover, {J'J^u)e{0) = e^^~^y^pej-^u so that 
/ := Jg^M satisfies /(O) G Cp, and therefore / G Hp(G). Note that by construction of 

the manifold, we have j-g^x ^ ~ i'e^(^)- 
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3.2. Closeness assumptions. Let us start this subsection with a lower bound on the 
operators H and in terms of the model parameters; for the definitions of the constants 
Ci/25 ^1/2 and Ci/2 see Lemma 12.21 and Lemma 12.101 Note that C1/2 still depends on 
ligiloo and 

Lemma 3.1. For e G (0,ei/2] the operators and H are hounded from below by Xq := 
— Ci/2. Moreover, if all lengths are finite, i.e. ie < 00, and q{v) < 0, then we have 

mfa{H)<^^^ and Ma{H,)< ^^""^ 



vol Xe vol Xe + £ vol X„ 

where Xe '■= Ue^^e ^'^ union of the edge neighbourhoods. 

Proof. We have to calculate the maximum of C1/2 and C1/2. Due to (13.21) we have 
|pp = voldXy and \q{v)\ = | Qdx„| < ||(5||oo volX„ so that 

Cy, < max{4cLllQ||L, < Cy, = max{64cLl|g||L, (3.7) 

where Cvoi := volX^,/ voldX^. The spectral estimates then follow by inserting suitable 
test functions into the Rayleigh quotients f)(/)/||/||^ and [)e(M)/||up. For /, we choose 
the edgewise constant function fe{x) = Pe- Note that / G H^{G). On the manifold, 

we choose the constant u := J^f = The upper bound on the infimum on the 

spectrum follows by the relation i^pl = volXg using (13.21) . □ 

Now we are in position to demonstrate that the two Hamiltonians are close to each 
other. We start with estimates of the identification operators and the forms f), f)e in 
terms of the "free" quadratic forms d and O^: 

Lemma 3.2. The identification operators J, J' = J* , , J'^ and the quadratic forms 
i)s and t) fulfil the estimates 

\\Jf - J'ff < SUfW^HG). \\J'^ - J"^f < ^sMhhx.), (3.8a) 

ii^/ir = ii/ip, \\j'ur<\\ur, (3.8b) 

J'J/ = /, \\JJ'u-uf<6'Ju\\f,,^j,^^, (3.8c) 

\[){J'\ f) - Uu, JV)| < <5.||M||HMx.)||/i|Hi(G) (3.8d) 
with 6s = 0(£:^/^) as e 0, being given explicitly by 

-2^1 / 2 



2 foeCvol £ A 2\ ^ , /1 , 



io 'A2(i?)' IX^iv) ' ' ioX2iv) 

2£/ 2 X 4£c.oi||g||L l (o 

Here, Cq := min{l,£_} = mine{l,4} < 1? ^2{.E) := miUe A2(e) andc^oi = volX^,/ voldXy. 
Moreover, A2(e) and X2{v) denote the second (first non-vanishing) eigenvalue of the 
(Neumann- )Laplacian on and X^, respectively. 

Proof. The first condition in (I3.8ap is here 

11^/ - J'fWk =eyo\XMv)\' < ^c.oi(||/1l| + fjlfWl) 

using Lemma \2A] with a = io and the fact that |pp = vol9X^ due to (13. 2p . Next we 
need to show the second estimate in (I3.8al) . In our situation, we have 

WJ'u - ru\\l = e'~' E WXeWlpllfu - fu{0)\' <e{l + ^) ||d«||^^,„ 
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using Lemma I^TSl with a = a = Eq. Moreover, (]3.8bp and the first equation in (I3.8cl) are 
easily seen to be fulfilled. The second estimate in fl3.8cp is more involved. Here, we have 

\\jj'u — uP=y IIm — -f -ull^ + IIm|| V • 

II II / J II J e ll^e.e II ll^e,i! 

e 

The first term can be estimated as in fl2.14l) by 

where u{s) := u{s, ■). The second term can be estimated by Lemma [2.91 so that 



-2 

0, > max< Ae 



1 / 2 \] £^ SeCvoH 
+ Cvol 1 H ' ^ 



which is sufficient for the estimate flS.ScI) . 

Let us finally prove fl3.8d|) in our model. Note that this estimate differs from the ones 
given in |P06] by the absence of the potential term = e~^Q there. In our situation, 
we have 

\i){j"u,f)-U^,j'f)\' 

< 2e'-^ - /e^(0)) (X'e, /')/.f + \q{v)ju - {Qu, l.)x„ri/(^)P 



Note that the derivative terms cancel on the edges due to the product structure of the 
metric and the fact that dy^-te = and the vertex contribution vanishes due to dx„l = 0. 
The first term can be estimated by 

2 \1 o 2e / 2 



2£ 5 + — -- - d« i < - 1 + 



aA2(t;) / a" "^^-^ " 4 V 4A2(^^)- 

using Cauchy-Schwarz, Lemma [2^ and the fact that HXellL ~ ^/^ — ^/^^ by our choice 
of Xe- For the second term, we use our definition q{v) = f-^ Q dx^ and q{v)fjl = 
{u,j^Ql^)x^ to conclude 



where P^u := u — f u is the projection onto the orthogonal complement of 1^,. The 
last estimate follows from fl2.14p . Collecting the error terms for the sesquilinear form 
estimate, we obtain 



1 



(5, > max< — 1 



2£A , 2 \ 4£c™i||Q|| 



2 

I oo 



4A2M/' 4A2M 



as lower bound on 5^, using also Lemma I^TT] for the estimate on and \\Q\\x ^ 

volX.IIQII^. " □ 

Now we can prove our main result on the approximation of a 5-coupling in the manifold 
model; for more details on the notion of "5-closeness" we refer to |P06| App.]. The 
resolvent estimate at 2; = i will be needed in Section H] when the lower bound Aq depends 
on e and may tend to — 00 as e — 0. Recall the definition of C1/2, < £1/2 (see (12.151) ) 
and Ao := — C1/2, and that C1/2 > C1/2. 
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Theorem 3.3. For e E (0,6:1/2]; ihe operators — Aq and H — \q are \plbe-close 
with 5e = 0(e^/^) given in (13. 9I) .- in other words, there is an identification operator 
J: l^iG) — > such that J* J = id, 

{id- J J*){H, - Ao)"M| < V26, and \\j{H - Xq)-^ - {H, - Ao)"V|| < 3^2(5,. 



Moreover, fore G (0,61/2] we have the estimate 

\\J{H -i)-^ -{H,-i)-'^j\\ < 104max{C'i/2, v^}, 
where ||-|| denotes the operator norm for operators from \-2{G) into L2{Xi;). 
Proof. The closeness of the operators H — Xq and H,, — Aq follows from the estimate 

\\du\\k + Mk < 2(f).H + (1 - Ao)||n||^J = 2\\iH, - Ao + lY^Mk 
by Corollary I2.1H and similarly for i7 on G by Corollary 12.31 and (13.71) , together with 
Lemma 13.21 The resolvent estimate can be seen as follows: Let R := {H — and 
R := {He - and let / G L2(G), u G l2iXe). Setting f ■= Rf e domH and 
u := RsU G domifj, we have 

{u,{JR-ReJ)f) = {u,Jf)-{u,Jf) 

= {u, (J - j')f) + {Uu, J'f) - W'u, /)) + {{J" - r)u, J) 

-i{{u,{j'-j)f) + {{j''-r)uj)), 

and therefore ^ _ _ 

\{u,iJR-R,J)f)\ < 105,max{Ci/2,v^}||/||||M|| 

using Lemmata 12.41 and I2.12[ and the fact that C1/2 < C1/2. □ 

Using the abstract results of |P06[ App. A] or |P08j . we can show the resolvent con- 
vergence and the convergence other functions of the operator: 

Theorem 3.4. We have 

\\J{H - - {He - z)-^J\\ = 0{e^/^), (3.10a) 

\\J{H - z)-^J' - {He - z)-^\\ = 0(5^/2) (3.10b) 

for z ^ [Aq, C)o). The error depends only on 6e, given in (13. 9p . and on z. Moreover, we 
can replace the function (p{X) = {X—z)~^ by any measurable, bounded function converging 
to a constant as X ^ 00 and being continuous in a neighbourhood of a{H). 

The following spectral convergence is also a consequence of the 0(e^/^)-closeness; for 
details of the uniform convergence of sets, i.e. the convergence in Hausdorff-distance 
sense we refer to |HN99i App. A] or jP08j . 

Theorem 3.5. The spectrum of He converges to the spectrum of H uniformly on any 
finite energy interval. The same is true for the essential spectrum. 

Proof. The spectral convergence is a direct consequence of the closeness, as it follows 
from the general theory developed in |P06l Appendix] and |P08] . □ 

For the discrete spectrum we have the following result: 

Theorem 3.6. For any X G crdisc(-f^) there exists a family {Xe}e with Xe G crdisc(-f^£) such 
that Xe ^ X as e 0. Moreover, the multiplicity is preserved. If X is a simple eigen- 
value with normalised eigenfunction (f, then there exists a family of simple normalised 
eigenf unctions {ife}e of He (e small) such that 

\\JlP - (fe\\x, 

as e ^ 0. 
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We remark that the convergence of higher-dimensional eigenspaces is also valid, how- 
ever, it requires some technicalities which we skip here. 

To summarise, we have shown that the 5-coupling with weighted entries can be ap- 
proximated by a geometric setting and a potential located on the vertex neighbourhood. 

Let us briefly sketch how to extend the above convergence results Theorems I3.3H3.6I 
to more complicated — even to non-compact — graphs. Denote by G a metric graph, 
given by the underlying discrete graph (V, E, d) with d: E — > V x V , de = ((9_e, d+e) 
denoting the initial and terminal vertex, and the length function i: E — > (0, oo), such 
that each edge e is identified with the interval /g = (0, ig) (for simplicity, we assume 
here that all length are finite, i.e., < oo). Let be the corresponding approximating 
manifold constructed from the building blocks e = h >^ and ,„ = eX^ as in 
Section |2.2[ For more details, we refer to |EP05l [P06l IEP081IP08] . Since a metric graph 
can be constructed from a number of star graphs with identified end points of the free 
ends, we can define global identification operators. We only have to assure that the 
global error we make is still uniformly bounded: 

Theorem 3.7. Assume thatG is a metric graph andX^ the corresponding approximating 
manifold constructed according to G. If 

inf X2{v) > 0, sup < oo, sup HQf v lU < oo, inf A2(e) > 0, inf 4 > 0, 

vev yizv vol oXy 

then the corresponding Hamiltonians H = Aq + '^yq{v)6tj and Hs = + J2v^~^Q^ 
are Ss-close, where the error Ss = 0{e^^^) depends only on the above indicated global 
constants. 



4. Approximation of the 5g-couPLiNGS 

The main aim of this section is to show how a the symmetrised (^'-coupling, or Si,, can 
be approximated using manifold model discussed above. To this aim we shall use a result 
of [CE04j by which a 5g-coupling can be approximated by means of several 5-couplings 
on the same metric graph, located close to the vertex and "lift" this approximation to 
the manifold. For the sake simplicity we will again consider the star-shape setting with 
a single vertex. We want to stress, however, that the method we use can be directly 
generalised to more complicated graphs but also, what is equally important, to other 
vertex couplings, once they can be approximated by combinations of (5-couplings on the 
graph, possibly with an addition of extra edges — see |ET06l [ET07j . 

Let thus G = be a star graph as in Section [2] where we denote the vertex in the 
centre by Vq and where we label the n = deg v edges by e = 1, . . . , n. Again for simplicity, 
we assume that all the (unsealed) transversal volumes pi = vol Yg are the same; without 
loss of generality we may put vol Yg = 1- Moreover, we assume that all lengths are finite, 
i.e. ie < oo, and equal, so we may put £e = 1. First we recall the definition of the 6'^- 
coupling: the operator H^, formally written as = Aq + PS'^^, acts as {H^ f)e = —f'l 
on each edge for functions / in the domain 

domi/^ := I / G HLxlG) I Vei,e2: 4(0) = 4(0) =: /'(O), ^/e(0) = /^/'(O), 

e 

Ve:/^(4)=0}. (4.1) 



For the sake of definiteness we imposed here Neumann conditions at the free ends of the 
edges, however, the choice is not substantial; we could use equally well Dirichlet or any 
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other boundary condition. The corresponding quadratic form is given as 

e e 

if /? ^ and 

[)^(/) = E ll/ell', domf)^ = { / G HL.(G) I E /e(0) = } 

e e 

if /3 = 0; the condition / G is obviously dual to the free (or Kirchhoff ) vertex coupling 
— see, e.g., p<u04l Sec. 3.2.3]. 

The (negative) spectrum of iJ^ is easily found: 

Lemma 4.1. If P > then > 0. On the other hand, if P < then has exactly 
one negative eigenvalue A = — where n is the solution of the equation 

(3k 

coshK + - sinhK = 0. (4-2) 

degf 

Proof. The non-negativity of if ^ follows from the quadratic form expression for /3 > 
and /3 = 0. We make the ansatz 

/e(s) = COshK(l — s) 

fulfilling automatically the Neumann condition at s = 1 and the continuity condition at 
s = since /e(0) = — KsinhK is independent of e. The remaining condition at zero leads 
to the above relation of k and /3, showing in another way that if /? > there cannot exist 
a negative eigenvalue. □ 

The main idea of the approximation of a (5s-coupling by Schrodinger operators on 
a manifold is to employ a combination of 5-couplings in an operator one may call an 
intermediate Hamiltonian H^'"', and then to use the approximations for (5-couplings given 
in the previous section. 

In order to define H^'"", we first modify the (discrete) structure of the graph G inserting 
additional vertices Vg of degree 2 on the edge e with the distance a G (0, 1) from the 
central vertex vq (see Figure [1]). Each edge e is splitted into two edges Ca and ci. 
We denote the metric graph with the additional vertices Ve and splitted edges by Ga, 
i.e., V{Ga) = {vq} U {ve\e = 1, . . . , n }, E{Ga) = { e^, ei | e = 1, . . . n } and £e„ = a, 
£ei = 1 — ct- This metrically equivalent graph Ga will be needed when associating the 
corresponding manifold. 

Remark 4.2. It is useful to note that the Laplacians and A^^ associated to the 
metric graphs G and Ga are unitarily equivalent. Indeed, introducing additional vertices 
of degree two does not change the original quadratic form dc with the domain H^(G) = 
doml) associated to the free operator A^^ = ii(G,o)- Figuratively speaking, the free 
operator does not see these vertices of degree two. We just have to change the coordinate 
on the edge e, i.e. we can either use the original coordinate s G (0, ie) on the edge e or 
we can spht the edge e into two edges Ca and ei of length iga = and £ei = — a = 1 — a 
with the corresponding coordinates. 

The core of the approximation lies in a suitable, a-dependent choice of the parameters 
of these ^-couplings. Writing the operator in terms of the formal notation introduced 
in (12. 6p . we put 

ii^''^ := A^ + h{a)5,, + V c(a)5,^, 6(a) = c{a) = --, 
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to be the intermediate Hamiltonian. Notice that the strength of central (5-couphng 
depends on j3 while the added ^-interactions are attractive, the sole parameter being the 
distance a. The operator can be defined via its quadratic form 

f)^'"(/) := E ll/eir - 4/(0) - - El/^WI'' domr = Hi(G), 

e e 

where H^(G) = Hp(G') with p = (1, . . . , 1), i.e. the functions / G H^(G) are distinguished 

by being continuous at vq, /ei(0) = /e2(0) =: /(O). 

The next theorem shows that the intermediate Hamiltonian converges indeed to the 
5g-coupling with the strength P on the star-shaped graph: 

Theorem 4.3 (Cheon, Exner). We have 

- z)-^ - {H" - z)-^\\ = 0(a) 

as a ^ Q for z ^M., where ||-|| denotes the operator norm on L2(G')0 

Note that the choice of the parameters h{a) c(a) of the 5-interactions as functions of 
the distance a follows from a careful analysis of the resolvents of H^'"" and if^. Each of 
these is highly singular as a — >■ 0, however, in the difference all the singularities cancel 
leaving us with a vanishing expression. Needless to say, that such a limiting process is 
highly non-generic. 




Figure 1. The intermediate graph picture used in the 5g-approximation 
and the corresponding manifold model. 

Let us now consider the manifold model approaching the intermediate situation Hamil- 
tonian H^'"^ in the limit £ — >■ with a = Oe = and < a < 1 to be specified later on. 
Let be a manifold model of the graph G as shown in Figure [H For the additional 
vertices of degree two we choose the vertex neighbourhoods as a part of the cylinder of 
length e and distance of order of from the central vertex vq. The edge =: now 
has the length = depending on e. The "free" edge ei joining with the free end 
point at s = 1 is again e-depending, namely it has the length 1 — = 1 — e". By the 
argument given in Lemma 12.71 we can deal with this error and assume that this edge 
again has length one, the price being an extra error of order 0(e"), affecting neither the 
final result nor the quantitative error estimate. Next we have to choose the potentials 



■^The claim made in |CE04| is only that the norm tends to zero, however, the rate with which it 
vanishes is obvious from the proof. We remove the superfluous degv from the definition of H^'"' in that 
paper. It should also be noted that the proof in |CE04| is given for star graphs with semi-infinite edges 
but the argument again modifies easily to the finite-length situation we consider for convenience here. 
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in the vicinity of the vertices v = vq and v = Vg- The simplest option is to assume that 
they are constant, 

e vol 

so that Qs^vdx = e~^qe{v) (see (12.101) and (13. 2p ). where we put 
q,{vo) := 6(£") = and := c(e") = 

The corresponding manifold Hamiltonian and the respective quadratic form are then 
given by 

= ^x. - ^"""^1^.0 - ^-'-^ E (4-3) 



and 



respectively. Note that the unsealed vertex neighbourhood X^^ of the added vertex 
has volume 1 by construction. 

Before proceeding to the approximation itself, let us first make some comments about 
the lower bounds of the operators if^'" and their manifold approximations iff: 

Lemma 4.4. If f3 < 0, then the spectrum of H^'"" is uniformly bounded from below as 
a —>■ 0, in other words, there is a constant C > such that 

inf a(ii^'") > -C as a ^ 0. 

-(f/? > 0, on the other hand, then the spectrum of H^'"" is asymptotically unbounded from 
below, 

inf (T(/i^''*) ^ -oo as a ^ 0. 

Note that although we know the limit spectrum as a ^ (see Lemma 14.11) . the re- 
solvent convergence of Theorem 14.31 does not necessarily imply the uniform boundedness 
from below of (see Remark 14. 8p . 

Proof. Let /3 < 0. Then an eigenf unction on the (original) edge e has the form 

y4cosh(Ks) + i?e sinh(Ks), < s < a 
CeCosh(K(l — s)), a < s < 1. 



fe{s) 



for K > 0, the corresponding eigenvalue being A = — k^. The Neumann condition /e(l) = 
at s = 1 is automatically fulfilled, as well as the continuity at s = for the different 
edges e, since /e(0) = A is independent of e. The continuity in s = a and the jump 
condition in the derivative lead to the system of equations 

A cosh(fi;a) + Be sinh(/ta) — Cg cosh(/t(l — a)) = 

— CeCosh/t(l — a) = k(— y4 sinh(Ka) — i?eCosh(Ka) — CeSinh/t(l — a)) 
a ^ ' 

e 

With the permutational invariance in mind, let us first analyse the situation with sym- 
metric coefficients. A, B = B^, C = Ce- Then X^e-^e ~ ^"^^ corresponding 
coefficient matrix for A, B and C vanishes iff 

-^(sinh(fi;a) coshK(l — a) — ok cosh + nniKasmhn — cosh(Ka) coshfi;(l — a)) = 
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leading to an eigenvalue A = — ^(a)^ of multiplicity one. It can be seen that K,{a) is 
bounded, and that the above equation reduces to (14.21) as a ^ oo. The other eigenvalues 
can be obtained from B and C as follows: set := e^'^'/". Then for k = 1, . . . ,n — 1, 
we have the coefficients -Be,fc = Q^^B and C(,,k = ©n'^C, e = 1, . . . , n. Since 

n n—1 

^i?,,, = 5^ef = 

e=l e=0 

for A; = 1, . . . , n — 1, we finally arrive at a coefficient matrix similar to the previous one, 
but with n replaced by zero. Consequently, if there were additional negative eigenvalues 
A = — K(a)^, they would be of multiplicity n — 1 and given by the relation 

sinh(/s:a) coshK(l — a) — Kacosh/t = 0. 

But this equation has no solutions for < a < 1 and k > 0. We skip the proof of this 
fact here. 

For the second part, assume that /5 > 0. It is sufficient to calculate the Rayleigh 
quotient for the constant test function / = 1 G H^(G) which yields 



being of order 0(a~^) if /? < and of order 0(a~^) if = 0, negative in both cases; 
recall that n = degf . □ 

Similarly, we expect the same behaviour for the operators on the manifold. 

Lemma 4.5. 1/(3 > 0, then the spectrum of is asymptotically unbounded from below, 
i.e., 

inf cr(iff ) — > — oo as e — >■ 0. 

Proof. Again, we plug the constant test function u = 1 into the Rayleigh quotient and 
obtain 

i)^{u) /3e-^^ + 

~ n{l + e + e°')+evo\Xyg 
which obviously tends to — oo as £ — 0. □ 

Remark 4.6. As for a counterpart to the other claim in Lemma 14. 4[ the proof of the 
uniform boundedness from below as e — for /? < seems to need quite subtle estimates 
to compare the effect of the two competing potentials on X^^vo and X^^^^ having strength 
proportional to and e"", respectively. Since the positive contribution Qe^vo = 

|^|^-i-2a more singular than the negative contributions Qe,ve = —£~^~", we expect 
that the threshold of the spectrum remains bounded as e — > 0. 

We can now prove our second main result. For the (5g-coupling Hamiltonian if^ and 
the approximating operator defined in f l4.ll) and fl4.3p . respectively, we make the 
following claim. 

Theorem 4.7. Assume that < a < 1/13, then 

\\{H^ -i)-^J - J{H^ ^0 

as e — > 0. 

Proof. Denote by H^'^ = H^'"-^ the e-depending intermediate Hamiltonian on the metric 
graph with (5-potentials of strength depending on e as defined before. For the corre- 
sponding graph and manifold model, the lower bound to lengths depends now on e. 
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specifically, £_ = = 5". Moreover, from the definition of the constants C1/2 < C1/2 
and £1/2 in (I2.15P and from Lemma [3.21 we conclude that 

= Cr,2{e) = 0(£-^"), £1/2 = £1/2(6) = 0(£3-) and 5 = 5, = 0{e^^-''^^/^). 

Note that the dominant term in the closeness-error 5, (see f l3.9p ) is the last one containing 
the potential. From Theorem 13.31 it follows now that 

||(iyf - i)-V - J(if^'^ - < 105, max{C7i/2(£), ^2} = 0{e^^-^^°'^'^). 

so that Theorem 14.31 yields the sought conclusion. Note that the exponent of e in 5eCii2{s) 
is (1 - 5a)/2 - 4a = (1 - 13a)/2 > provided < a < 1/13. □ 

We can now proceed and state similar results as in Theorems I3.4H3.7I for the 5'^- 
approximation by using arguments similar to those in [P06| App.] or [P08] . where only 
non-negative operators were considered (covering, as usual, operators bounded uniformly 
from below by a suitable shift). In our present situation, we can only guarantee the 
resolvent convergence at non-real points like z = i. Nevertheless, the arguments in [P06t 
App.] or |P08j can be used to conclude the convergence of suitable functions of operators 
as well as the convergence of the dimension of spectral projections. 

Remark 4.8. Note that the asymptotic lower unboundedness of (and of the interme- 
diate operator H^''^) for /3 > described in Lemmata 14.41 and 14.51 is not a contradiction 
to the fact that the limit operator is non-negative. For example, the spectral conver- 
gence of Theorem [23] holds only for compact intervals / C R. In particular, a{H^)r\I = 
implies that 

a{H^) n / = and a{H^'') H / = 

provided £ > is sufficiently small. This spectral convergence means that the negative 
spectral branches of all have to tend to — cx3. 

Acknowledgement. O.P. enjoyed the hospitality in the Doppler Institute where a part 
of the work was done. The research was supported by the Czech Ministry of Education, 
Youth and Sports within the project LC06002. 

References 

[ACF07] S. Albeverio, C. Cacciapuoti, and D. Finco, Coupling in the singular limit of thin quantum 

waveguides, J. Math. Phys. 48 (2007), 032103. 
[CE07] C. Cacciapuoti and P. Exner, Nontrivial edge coupling from a Dirichlet network squeezing: 

the case of a bent waveguide, J. Phys. A 40 (2007), L511-L523. 
[CE04] T. Cheon and P. Exner, An approximation to 5' couplings on graphs, J. Phys. A 37 (2004), 

no. 29, L329-L335. 

[CS98] T. Chcon and T. Shigehara, Realizing discontinuous wave functions with renormalized short- 
range potentials, Phys. Lett. A 243 (1998), 111-116. 

[DEKOl] P. Duclos, P. Exner, and D. Krejcifik, Bound states in curved quantum layers. Comm. Math. 
Phys. 223 (2001), no. 1, 13-28. 

[DT06] G. F. Deir Antonio and L. Tenuta, Quantum graphs as holonomic constraints, J. Math. Phys. 
47 (2006), no. 7, 072102, 21. 

[E97] P. Exner, A duality between Schrddinger operators on graphs and certain Jacobi matrices, 

Ann. Inst. H. Poincare Phys. Theor. 66 (1997), no. 4, 359-371. 

[EKK+08] P. Exner, J. P. Keating, P. Kuchment, T. Sunada, and A. Teplayaev (eds.). Analysis on 
graphs and its applications, Proc. Symp. Pure Math., voL 77, Providence, R.I., Amer. Math. 
Soc, 2008. 

[ENZOl] P. Exner, H. Neidhardt, and V.A. Zagrebnov, Potential approximations to 5' : an inverse 
Klauder phenomenon with norm-resolvent convergence, Comm. Math. Phys. 224 (2001), 
no. 3, 593-612. 

[EP05] P. Exner and O. Post, Convergence of spectra of graph-like thin manifolds. Journal of Ge- 
ometry and Physics 54 (2005), 77-115. 



APPROXIMATION OF VERTEX COUPLINGS BY SCHRODINGER OPERATORS ON MANIF0LDS19 



[EP07] , Convergence of resonances on thin branched quantum wave guides, J. Math. Phys. 

48 (2007), 092104+43. 

[EPOS] , Quantum networks modelled by graphs, Quantum Few-Body Systems, AIP Conf. 

Proc, vol. 998, Amer. Inst. Phys., Melville, NY, 2008, pp. 1-17. 

[ET06] P. Exner and O. Turek, Approximations of permutation-symmetric vertex couplings in quan- 
tum graphs. Quantum graphs and their applications, Contemp. Math., vol. 415, Amer. Math. 
Soc, Providence, Rl, 2006, pp. 109-120. 

[ET07] , Approximations of singular vertex couplings in quantum graphs, Rev. Math. Phys. 

19 (2007), no. 6, 571-606. 

[FW93] M. I. Freidlin and A. D. Wentzell, Diffusion processes on graphs and the averaging principle, 
Ann. Probab. 21 (1993), no. 4, 2215-2245. 

[G08] D. Gricscr, Spectra of graph neighborhoods and scattering, Proc. London Math. Soc. (3) 97 
(2008), no. 3, 718 752. 

[HN99] I. Herbst and S. Nakamura, Schrodinger operators with strong magnetic fields: Quasi- 
periodicity of spectral orbits and topology, American Mathematical Society. Transl., Ser. 2, 
Am. Math. Soc. 189(41) (1999), 105-123. 

[Ku04] P. Kuchment, Quantum graphs: I. Some basic structures. Waves Random Media 14 (2004), 
S107 S128. 

[KuZOl] P. Kuchment and H. Zcng, Convergence of spectra of mesoscopic systems collapsing onto a 

graph, J. Math. Anal. Appl. 258 (2001), no. 2, 671 700. 
[MV07] S. Molchanov and B. Vainbcrg, Scattering solutions in networks of thin fibers: small diameter 

asymptotics, Comm. Math. Phys. 273 (2007), no. 2, 533-559. 
[Pa07a] B. Pavlov, Neumann Schrodinger 2D junction: collapsing on a quantum graph: a generalized 

Kirchhoff boundary condition. Preprint (2007) . 
[Pa07b] , A star-graph model via operator extension, Math. Proc. Cambridge Philos. Soc. 142 

(2007), 365-384. 

[P05] O. Post, Branched quantum wave guides with Dirichlet boundary conditions: the decoupling 

case. Journal of Physics A: Mathematical and General 38 (2005), no. 22, 4917-4931. 

[P06] , Spectral convergence of quasi- one- dimensional spaces, Ann. Henri Poincare 7 (2006), 

no. 5, 933-973. 

[P08] , Spectral analysis on graph-like spaces, Habilitation thesis in preparation (2008). 

[RSOl] J. Rubinstein and M. Schatzman, Variational problems on multiply connected thin strips. I. 

Basic estimates and convergence of the Laplacian spectrum. Arch. Ration. Mech. Anal. 160 

(2001), no. 4, 271-308. 

[SaiOO] Y Saito, The limiting equation for Neumann Laplacians on shrinking domains.. Electron. J. 
Differ. Equ. 31 (2000), 25 p. 

Department of Theoretical Physics, NPI, Academy of Sciences, 25068 Rez near Prague, 
AND DOPPLER Institute, Czech Technical University, Brehova 7, 11519 Prague, Czechia 
E-mail address: exner@ujf.cas.cz 

Institut fur Mathematik, Humboldt-Universitat zu Berlin, Rudower Chaussee 25, 
12489 Berlin, Germany 

E-mail address: post@math.hu-berlin.de 



